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' It is demonstrated how a convenient choice of the mathematical 

^ . structure of the quantum cosmology superspace, precisely the defini- 

tion of a convenient regular state superspace and the restriction of the 
dynamics to this space, yields directly to an irreversible evolution, in 
the classical (and semiclassical) phase of the universe, where: 



U 

• Decoherence and correlations take place and therefore give origin 
to a classical universe. 

^ ■ • The second law of thermodynamic is demonstrated. 

' • Connection with Reichenbach branch system idea can be imple- 

mented. 

Some rough coincidence with observational data are obtained. 
The arrows of time can be correlated. 

Time asymmetry can be explained as a state space asymmetry ( 
e. g. like a spontaneous symmetry breaking) 
All these facts solve the problem of time-asymmetry and show 
that it is time asymmetry itself that defines the most important 
features of mathematical structure of superspace. 
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1 Introduction 



The role of physics is to explain nature in the best possible way. Therefore 
physicists use to consider a set of physical phenomena and to choose the 
best axiomatic structure to mimic these phenomena. This axiomatic struc- 
ture contains a mathematical structure and a set of axioms (or postulates, 
or principles, or hypotheses), stated using the language of the chosen mathe- 
matical structure. If the mathematical structure is the most naturally related 
to the set of phenomena and a minimal number of axioms is used, physicist 
say that they have explained nature in the best possible way. But, frequently 
a better mathematical structure and a smaller set of axioms are found to 
explain a larger set of physical phenomena, then physicists say that they 
understand the problem even better because, in fact, they have found a bet- 
ter explanation i. e. a more economical one. But the chosen mathematical 
structure and the chosen axioms cannot be explained by themselves, since 
the only motivation of the choice is to mimic nature in the best possible way. 

Thus, physical phenomena are not a consequence of the chosen math- 
ematical structure, quite on the contrary, the choice of the mathematical 
structure is a consequence of the physical phenomena that we are trying to 
explain. 

Gravitation was explained by Einstein choosing a riemannian manifold as 
mathematical structure and postulating that particles space-time paths were 
the geodesic of such a manifold, etc. The only explanation of this choice is 
that the theoretical motions, so described, mimic real motions better than 
the motions described using other curves or other kind of manifolds (flat 
space-time, projective manifolds, etc.). But the riemannian manifold and 
the geodesies cannot be explained by themselves. In fact, the choice of the 
riemannian manifold, as the mathematical structure to explain gravity, is 
really motivated by the method we use to measure time and distance ||T]|. 

It also happens that, when a new or unfamiliar mathematical structure 
is introduced, some physicists think that the new structure is introduced by 
hand, because they do not realize that every mathematical structure was 
introduced by hand in, order to explain nature in the best possible way. 

In the perspective of this pedagogical (perhaps pedantic but necessary) 
introduction we will choose: 
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1.1 The set of phenomena. 

The set of phenomena considered in this paper will be those of usual quantum 
cosmology (QC) based in Wheeler-De Witt (WDW) equation [|[, [|], plus 
those of irreversible statistical quantum physics, like the definitions of the 
various arrows of time, the final equilibrium state, decoherence, correlations, 
etc. (see the book |^ where almost all subjects, results, and problems can 
be found). 

1.2 The mathematical structure. 

QC is based in WDW equation; 

H^ = (1) 

where H is the hamiltonian operator of the model for the considered universe 
and \1/ the wave function of the universe (that will be also called 1^1/ >) 0], 
p. Usually H is well defined by the model we are studying, but the real 
mathematical nature of \1/ is not so well defined. In fact, let us call "g" the 
configuration variables and Q the configuration space, then \1/ = \l/(g) is a 
function of the configuration variables. We must define the space of wave 
functions Let us suppose that this space has a discrete basis V of certain 
functions of the q (let us say polynomials of the q or polynomials multiplied 
by convenient dumping factors for those q, s that goes to infinity). But, this 
set is not enough because usually wave functions belong to a complete space. 
Then we must complete the space expanded by V with some topology. The 
usual idea is to complete the space expanded by V with the topology of 
the norm. Even if the choice of a norm contains a very important physical 
problem: the definition of probability, we are now exclusively interested in 
the mathematical problem of how to complete V with an adequate topology. 
Then we have really two problems: 

i.- We must define a norm in the space expanded by V something like [^]: 

II ^ 11= / <i/*{q)-qf{q)dq (2) 
JQ 

But we can also foliate Q with hypersurfaces S and define 0: 

II ^ 11= j^^*{q)dn^{q)dan (3) 
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where n is the normal to S and dan is the hypervolume of the hypersurface 
E differential element (this quantity will be positive definite in a convenient, 
subspace of superspace) . Or we can consider that really is an operator and 
then we must go to a third quantization 0, etc. So we have many possibilities 
to define the norm. Actually we are trying to copy usual quantum mechanics 
and the essential property to define an adequate norm is that it would be 
a constant under time evolution. But in QC there is no time [§] [|10|] . 
In fact, eq. (|ID is not a Schroedinger equation but an eigen equation that 
defines a "stationary" eigenfunction, with no time variable in it. Therefore, 
as the definition of norm is intimately intertwined with the problem of the 
definition of time in QC, and this problem is not solved P], [|12] it is 



not possible to give a definitive definition of norm. Really in this paper we 
will adopt the conservative attitude that this problem cannot be solved and 
therefore QC is a timeless theory , . As a consequence it is very difficult 
to define a satisfactory norm. 

ii.-But even if a norm would be chosen, say (^, we have another problem. 
If we complete the space expanded by V with we will find H =L'^{Q) the 
Hilbert space of square integrable functions over Q. But H is an operator 
with derivatives and the functions of Ti are, in general, not derivable, they 
are just square integrable, so for an arbitrary function \1/ of 7^, has no 
meaning (and we would have the same problem with the other norms). So 
we must choose another topology to complete the space expanded by V. 
Let us consider the Schwarz class functions, namely functions that can be 
derived an infinite number of times and that are well behaved in the eventual 
infinite of coordinates q (precisely they vanish faster than any polynomial). 
The set of Schwarz functions is a nuclear space. Then if we complete the 
space expanded by V with the corresponding topology, which is not a norm 
topology but a nuclear one |T^, we obtain a space S where we can derive 
any number of times. If is the space of (anti)linear functional over 5, 
our mathematical structure is really the Gel'fand triplet (or rigged Hilbert 
space) : 

5 C 7^ C 5^ (4) 

Then we can say that S is the superspace of regular states where we can do 
all our computations and is the superspace of generalized states, where 
we will find usual distributions like Dirac's deltas or plane (or curved) waves, 
that can be used to expand regular states. In ordinary quantum mechanics 
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H would be the superspace of states, but in QC this space looses almost all 
its importance because we do not know what norm we have to use for sure 
(even if, to compute probabilities , we will use norm in spite of the fact 
that we have other possible choices [p!^] ), while S and x are well defined, since 
their definitions are norm-independent. The "inner product" < F|\E' > of an 
element of written as a "bra", < F\, with an element of S written like a 
"ket", >, is also well defined, since it is just the functional acting over a 
regular state, (also the product in the inverted order is well defined, if we set 
< \E'|-F >=< -F|\I' >*). So S is the arena where timeless QC works: where 
time and norm have disappeared from the mathematical structure and are 
substituted by the nuclear topology of S. 

But if we use 5 as a regular space we cannot encompass (in the most 
economical way) irreversible statistical quantum mechanics, which has an 
asymmetry that it is not contained in S. There are only two causes for 
asymmetry in nature: either the laws of nature are asymmetric or the so- 
lutions of the equation of the theory are asymmetric. E. g.: The laws of 
nature are asymmetric in the case of the weak interaction. The solutions of 
the theory are asymmetric in the case of spontaneous symmetry breaking. 

Time asymmetry is not an exception. Thus, if we want to retain the 
time-symmetric laws of nature (namely the symmetry of eq. (0)) the only 
way we have to explain the time-asymmetry of the universe or its subsystems 
is to postulate that the space of solutions is not time-symmetric, namely we 
use the second cause of asymmetry. So the proper way to solve the problem 
is simply to define a realistic time-asymmetric space of physical admissible 
solutions 0_, i. e. an adequate mathematical structure for superspace. 0_ 
will contain the states that evolve in an admissible way (e.g.: Gibbs ink 
drop spreading in a glass of water, a sugar lump solving in a cup of coffee, 
etc.) and will not contain the non admissible evolutions ( the ink or the 
sugar concentrating spontaneously and creating the drop or the lump). The 
problem is to choose in the best possible way. As we will see we will choose 
one with the required asymmetry. We will follow a heuristic approach: we 
will suppose that H is endowed with all the properties necessary to define 
a reasonable universe. Obviously these "realistic hamiltonians" are the only 
H that we must consider. Therefore practically, in H there must be always 
some fields, like the mater field, the electromagnetical field, and also the 
gravitational field (namely only the graviton field). So we would write this 
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hamiltonian in a midi-superspace way as: 



H = hnfigj, Txj) + hf{Lp,p^) + hi{gij^ip) (5) 

where hnf{gj,iTj) is the "no field" hamiltonian (a function of a discrete set of 
modes of the gravitational field gj and the corresponding momenta tcj), hf 
is the "fields" hamiltonian, (let us say the continuous set of all the modes of 
just one field (p , that represents all the physical fields in our model, and the 
corresponding momenta p^), and hi is the interaction hamiltonian among 
fields and no fields (usually an interaction among only the configuration 
variables) . 

To complete the panorama, (using the Gel'fand-Maurin theorem |T^) we 
can find, for the partial hamiltonian hnf + hf, a spectral decomposition: 

roc 

hnf + hf = "^uiili >< i\ + / >< Lj\duj (6) 

where the first term of the r.h.s. corresponds to the spectral decomposition 
of hnf, (cf. eq. ||) and therefore to a discrete spectrum, and the second one 
to the spectral decomposition of hf and therefore to a continuous spectrum. 
Some negative eigenvalues "u" must appear because H is not bounded from 
below. Any state |^ >G 5 can be expanded as: 

1^ >= ^ |z >< z|^ > + / \iu >< uj\'^ > du (7) 
. Jo 

This is the spectral decomposition of |\E' > in basis {|i >, |c<j >}. In order 
that > would satisfy the WDW eq. (|1|) it must be: 



roo 

>= VcjjU >< z|^ > + / uj\uj >< uj\'^ > duj + hi\-^ >-- 
A JO 







Most likely this equation can be solved. 

Also, in every reasonable model, the universe will end in a classical phase, 
going first through semiclassical one. Namely, a variable a will exist (one 
of the gj or a function of the gj ) such that a time t] can be defined as 
a function of a. When ^ oo we will obtain a classical geometry g"^^ 
for the universe. It will be the most probable geometry of the universe 
(i. e. the geometry that appears most frequently) [0. Using time t] we 



can transform eq. (|1|) in a Schroedinger equation, with the corresponding 
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hamiltonian h Then using h and the classical geometry g"^ we can find 
a semi-classical vacuum state |0,OMt > for the fields (the so called adiabatic 
vacuum), that diagonalizes the hamiltonian h = h{g^) = h{out) (computed 
in the geometry g°^^), the creation and annihilation operators related to this 
vacuum, and the corresponding Fock space. But the only essential ingredient 
we need to implement the theory is the h{out). Then, using these objects, 
we can find a set of eigenvectors \uj, out > [O], such that: 



h\uj,out >= uj\u!,out > (9) 
where a; is a continuous eigenvalue of h (say < u; < oo) So: 

roo 

h= uj\uj, out >< uj, ou\du (10) 
Jo 

where \uj,out >G , thus if |\1' >G S, < u!,out\'^ > is well defined and 
so /i|\E' >. The existence of this kind of expansions can be also considered a 
consequence of Gel'fand-Maurin theorem [|T^]. But all these manipulations 
are just formal, so we will be sure that what we are doing is correct only in 
concrete examples, as the one in the next section. In fact, we will find all 
this mathematical elements in the model presented there. Surely we will also 
find these elements in more complex models |jl9l . 



Now, let us define our new regular state space 0_ C S. Precisely, we 
can promote to a complex variable z and ask, not only that >G S, but 
also that < z,0Mt|\E' > would be an analytic function in the lower complex 
halfplane (precisely that < co',OMt|\E' >G being the Hardy functions 
class from below). These functions |\E' > belong to a space 0_ such that: 

0_ c 5 (11) 

and we have a new Gel'fand triplet: 

0_ C Ti: C 0!! (12) 



then we know that 13 



5^ c (13) 

So we have restricted the regular state superspace and simultaneously we 
have enlarged the generalized state superspace so we will have more general 
spectral expansions (this fact will be of outmost importance). But, as we 
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can as well choose the upper complex halfplane (precisely < cu, OMt|\E' >G 
being the Hardy functions class from above), we also have another space 
0+ such that: 

0+ C 5 (14) 

and another Gel'fand triplet: 

(j)+cnc(t)i (15) 

We also know that HTSi: 



5x C 0^ (16) 

Now we can also say that we have obtained the space 0_ completing the space 
expanded by V with the nuclear topology of 0_, namely the 5-topology re- 
stricted to 0_ (and the same thing can be said about 0+). Clearly this topol- 
ogy is endowed with a new asymmetry which S does not have. Precisely, 
this asymmetry allows us to choose between 0_ or 0+ even if we maintain 
all the symmetries of H [^. Thus we can break one of these symmetries, 
restricting the dynamics to superspace 0_ which then would be considered 
as the superspace of regular states. As we will see this restriction produces 
the desired time asymmetry. Frequently physicists make analytic continua- 
tion in the complex energy plane supposing that some functions are analytic 
in one halfplane only. In these cases they are implicitly using the kind of 
mathematical structure we have explicitly introduced here, so the idea is, by 
no means, new. 

Thus our mathematical structure will essentially be (|12D, 0_ will be our 
superspace of regular states, where we must find the states that satisfy the 
WDW eq. (|^) and 0^ will be our generalized state superspace. From (pT]) we 
see that we have restricted our regular state superspace, so nothing unphys- 
ical can happen. We are just adding a new requirement to regular states, 
in order to assure their asymmetry. Furthermore 0_ is dense in Ti, so if 
someone would believe that the real "physical" states belong to 7i we could 
argue that these "physical" states can be approximated, by regular states of 
0_, as close as we want. 

1.3 The axiomatic structure. 

We do not pretend to give a completely rigorous axiomatic structure in this 
paper (but just an approximation of it). Furthermore we do not know if the 



8 



proposed axiomatic structure is unique. We are just proposing a first draft of 
a complete axiomatic structure, so we will call to our axioms just hypotheses. 

As we would like that the equilibrium state would be contained in our 
theory we must also consider mixed states p and, therefore, the spaces: 

$_ = 0_ ® 0_, = n®n, = 0^0 0^ (17) 

and work in the Liouville triplet: 

C £ C $!! (18) 

where C is the usual Liouville space of ordinary mixed states, which actually 
we will never use, since our regular state superspace is So, let p be a 
selfadjoint density matrix. Then our main hypotheses are.- 
Hi. The state p of the universe satisfies the equations P^j: 

Hp = (19) 



H2. The state p of the universe belongs to the superspace i. e..' 

pe^^ (20) 



H^. p{q, q') is proportional to the correlation between the configurations q 
and q' and p{q, q) is proportional to the probability to find the configuration 
q in the universe. 

These three axioms correspond to the three elements necessary to go into 
QC mentioned in the introduction of paper 0: dynamics, "initial condition" 
(precisely definition of the physically admissible states of the universe in a 
timeless formalism), and interpretation. Of course H2 alone does not fix 
the actual state of the universe, but any state of the universe we choose to 
build our theory must be contained in if we want that time-asymmetry 
would appear in a natural way. In this paper we do not address the prob- 
lem to find the real and unique state of the universe, but only to define a 
superspace of admissible states such that the universe would turns out to be 
time- asymmetric. 

We will see how far we can go with this axiomatic structure. The paper 
is organized as follows: 
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• In section 2 we introduce our model, its semiclassical approximation 
and we obtain a new spectral decomposition, using the regular super- 
space of hypothesis H2. 

• In section 3 we obtain the evolution equation of the states. 

• In section 4 up to section 9 we find the physical characteristics of the 
model. 

• In section 10 we draw our main conclusions. 

• Two appendices, containing further explanations and the way to fix the 
two ± signs of the theory, complete the paper. 



2 The Model. 

Let us see how we can implement all we have said in a simple model. 

Let us consider the model of section 3, of ref. |22|], or better the one of 



ref. p3[, where a Robert son- Walker metric: 

ds^ = a\i]){dr]^ - dx^ - dy^ - dz'^) (21) 

is studied, (we will mostly consider the flat space geometry case). The total 
action is: S = Sg + Sf + Si, being Sg{= Snf) the gravitational action, 5"/ the 
usual action of an spinless massless field ip, conformally (^ = |) coupled , 
and Si the interaction given by a mass term in Robertson- Walker geometry. 
The gravitational action is given by: 

Sg = M' Jdr][~a -V{a)] (22) 

where M is the Planck's mass, r] is the conformal time, a is the Robertson- 
Walker scale , a= and V{a) is a potential that arises from the spatial 
curvature, a possible cosmological constant, and, eventually a classical matter 
field. V{a) is a potential with a bounded support contained in < a < Oi, 
with ai 3> 0, (in many examples V{a) is a function of and V{a) strongly 
vanishes, for ^ 00, so our potential can be considered as a good 



approximation of these examples). This case is the simplest of all, but we 
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believe that the main features, that we will find, will also be present in more 
general cases. 

The WDW eq. (|T]) for our model is: 

H^{a,ip) = {hg + hf + h,)^{a,^) = (23) 

where (in our the fiat space geometry case): 

K = l^A + MV(a) (24) 



h 



f 



\ lidl - kV,)dk (25) 



hi = ^J^ ^'rfk (26) 

where m is the mass of the scalar field, and k"^ = Ikl^, where - is the linear 
momentum of the field, in the fiat case we are working (in the two other cases, 
namely open and close space geometry, the integrals of eqs. (^) and (^) 
are integrations on adapted coordinates, in the open case, and summatories, 
in the close case, where k is substituted by a discrete variable; see the 



corresponding equations in [22 



Now, let us go to the semiclassical case using the WKB method [^, 
[g. So let: 

■^{a, ^) = exp[iMSia)]xia, ip) (27) 

To satisfy WDW eq. (|I|), at order M^, the principal Jacobi function S{a) 
must satisfy the Hamilton- Jacobi equation: 



'dsy 

da I 



2V{a) (28) 



Now we can define the time, in our up to now, timeless theory. It is the 
(semi) classical time parameter rj = ri[a) given by: 

d dS d _^ /2V( ) ^ (29) 
dr] da da ^ da 

From eqs. (pSf ) and (^) we can find the set of classical solutions: 

a = ±f{v,C) (30) 



11 



where C is an arbitrary integration constant. Using different values for this 
constant and different choices of the ± sign we obtain different classical 
geometries (in more general cases many constants would be necessary). For 

a > ai_ it is y2V(a) = (since V{a) have a bounded support, contained in 
[0,ai]), and we cannot define time r] using eq. (pOf), thus we must choose 
another hand for our clock to define the time there. To avoid this problem 



let us consider that when a > Oi it is y2V{a) = e = const. = 0. We can 
always make e = to reobtain. the real case. Then a will be 

a = ±eri + C (31) 

So we can see that the potential can also be considered as a function with 
bounded support in variable?7. We will consider that e > 0. The role of C 
is just to fix the origin of time, so we can take any C we want . As the 



coupling is conformal we will have well defined vacua , ||26[ . In particular 
we can consider two scales a^n and Qout such that < aj„ <C ai, aout ^ 
and define the \0,in >, \0,out > vacua there. (We can as well transform all 
the equations to the non-rescaled case, consider the proper time t = J adr], 
and the physical momentum - and define the |0, out > in the r] ^ oo limit, 



as in the appendix A of paper |2^, but here we will use the first simpler 
formalism.) 

In this model we have: 

h{a) = hf{ipk) + hi{a,ipk) (32) 
where we have omitted the in h{a). Then: 

= + ^^^^ 

where (cf. eqs. ( P^D and ([2^)): 

nl{a)=m'a^ + e (34) 

So h{a) is a time dependent hamiltonian, where all its time dependence comes 
from a scale variable mass rn^a?. It is well known (|2^, I^Bj) that we can 
diagonalize this time dependent hamiltonian at aj„ and at aout and define the 
corresponding vacua, the corresponding creation and annihilation operator. 
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and the corresponding Fock spaces. For the out-geometry the vacuum will 
be the adiabatic vacuum, since aout ^ oi, therefore all the out elements will 
coincide with those defined in the introduction. In fact, the out-geometry is 
almost constant during the final time evolution (that goes up to ^ oo) and 
therefore they correspond to the geometry with the maximum probability. 
h{aout) reads: 

Kaout) = / ^kia'out)aLt,k(^out,kdk (35) 

J k ' 

where a]^^ and aout,k are the creation and annihilation operator correspond- 
ing to the out vacuum. With these objects we can construct the Fock space 
with a basis: 

\ki.k2,...kn,out >= \{k},out >~ aL,fci«L,fc2-"«L,fcjO,OMt > (36) 

where we have called {k} the set ki, k2, ...kn- These states are eigenvectors 
of h{aout), precisely; 

h{aout)\{k},out >= u{aout)\{k},out > (37) 

(of course, this equation corresponds to eq. (^) where: 

^(aout) = J2 ^k{out) (38) 

fcG{fc} 

We can use this energy to label the eigenvectors as: 

\{k},out>=\uj,[k],out> (39) 

where [k] is the remaining set of labels necessary to define the vector unam- 
biguously. [k],out >} is an orthonormal basis, namely: 

< u, [k], out\u'[k'],out >= 6{uj - uj')6{[k] - [k']) (40) 

/•OO r 

1= / duj \uj,[k],out >< uj,\k],out\d[k] (41) 

JO J[k] 

where the meaning of the symbols related with [k] is evident. In the close 
space geometry case the indices would be discrete and the integral a summa- 
tory. 
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All the same can also be done at aj„. We can now define the S-matrix 
between the in and out states: 



\k], 



jfc/] =< u, [k], in\uj', [k'], out > 



(42) 



According to papers |2j, p8[, this matrix has an infinite set of complex poles 
as we will demonstrate in section 4 (also an example is given in |23| and using 
this paper and |^ and other examples can be obtained). 

If we forget the indices [k] and consider again the eq. (H), we see that 
\uj, [k],out > is the \uj,out > of this equation. In the introduction we have 
defined the triplets (|T^ and ([l5|) only using the hamiltonian h{out) = h{aout)- 
These triplets correspond to a Fock space defined for aout- But there will also 
be two similar triplets defined in the Fock space at a^^. We make the following 
choice (motivated on reasons that will be apparent in a moment): for the 
in- Fock space we will use functions Iv^ >G 0+m, namely such that \(p S 
and < u!,in\ip >G H"^, and for the out-Fock space we will use functions such 



ip >G S and < u!,out\ip >G H^. So 
and our regular states belong to 



,out 



that \(f >G (p-^out, namely such that 
the 0_ of the introduction is now 0. 
this space [^. The role of 0+,in is to allow us to define the corresponding 
functional space 0+ j„ so we can use the functionals of this space in some 
spectral decompositions PT|. As both vacua, at ai„, and aout are well defined 



and the particle production between these vacua is finite and the theory is 
implementable We can then multiply the state of both Fock spaces. 
So let us again write eq. (0) with no [k] : 



duj\uj, out >< uj, out\ 



(43) 



Of course there is an analogous equation for the "in" case. Now using this 
equation and eq. (H^) we have: 



roc roo 

\uj,out>= / duj'\uj' ,in >< uj' .in\uj, out >= / duj'\uj',in > S^' ^ (44) 
Jo Jo 



Then: 



oo POO 



duj I duj\uj\in > S^^i < iJ,out\ 

JO 



or: 



oo poo 



(45) 



< >= / duJ du' < ip\u' , in > S^^' ^ < u , out\(p > (46) 
JO Jo 
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Now let Iv^ >G 4>-,out and \ip >G 0+,m- Then < z, out\ip >G H^, < z, inlip >G 
and therefore < ip\z,in >G H^. So in the integrand of the last equation 
all the factors are analytic in the lower halfplane, with the exception of 5'(^,aj', 
which has an infinite number of poles Qn as we have already said. Then we 
can choose any curve F, beginning at the origin, and going below all the poles 



of the lower halfplane up to the infinity of the positive real axis |^3|. We can 
now change the integration contour of eq. (|46D from [0, oo) to the curve F . 
If we add the pole contribution as in papers [0, p5[, and p6|], we obtain: 



< iplip >-- 



E 

n 



< tp\n+ >< n—\ip > + J dz J dz < tjj\z',in > Sz',z < z,out\(f >- 



V < ilj\n+ >< n - \ip > + I dz < %Ij\z+ X z-\<^> 



(47) 



where the summatory comes from the residues of the poles (each pole Qn is 
labelled by a discrete index n of course ImQn < 0) and: 



\z, — >= \z,out >, \z+ >= J dz' Sz'z\z' ,in > (48) 

Then, in a weak sense, we have found a new spectral decomposition of 1: 

1 = E X ^ ^ I + / dz\z+ X z - I (49) 

Following the same procedure with < ip\h{out)\(f > we can obtain the spec- 
tral decomposition of h{out) (always in a weak sense): 



h{out) = >< n — \+ z\z+ >< z — \dz 

We have three possibilities to choose the curve F: 



(50) 



i.- To use all possible curves F as in paper [36 



ii. - To take the curve (-oo, 0], in the second sheet, as in ref. provided 
we have a good behavior at infinity of the lower halfplane. 

iii. - To use the Nakanishi trick |3^, as in paper ||3^, namely to define 
tilded functionals such that: 



\z+ >< z — \dz 



\uj+ >< U!—\dtO 



(51) 
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We will use this last notation. Then we have: 



1 = V \n+ >< n-\+ duj\uj+ >< to - \ (52) 

/i(o-ut) = >< n — I + / uj\uj+ >< uj—\du; (53) 

„ JO 



(see also |]38[ and paper [Q, where the Nakanishi trick is explain). From its 
own definition it is evident that \n— >, \z— >, \uj— >G 
vectors are functionals over \(f >G 



^ out, since these 



,out 



and that \n+ >, \z+ >, \uj+ >e 
(f)^^in since these vectors are functionals over {ip >G Restoring the [k] 

and eliminating the tilde (as we will always do below) the last equation reads: 

h{out) = V fi„|n+ >< n - I + / duo l VLk\uo, [k]+ >< uj, [k] - \d[k] (54) 
„ Jo Jk 

but we will continue with the previous shorthand notation and we will not 
write the [k] any more. It can be proved that the bases {|^+ >, \oJ+ > 
}, {\n— >, \uj— >} are a biorthonormal system, ||36|, and ||3^ namely: 



< n — \n'+ >= 6n 



< UJ 



< n 



\uj- 



>= 



(55) 



(56) 



\n+ >= 0, ... < UJ - \uj'+ >= 6{uj - uj') 
From all these equation we have that: 

h{out)\n+ >= Qn\n+ > 

< n — \h{out) = Vtn < n ~ \ 

where f2„ is a complex eigen value, and |n+ > are right eigenvectors, and 
<n—\ left eigenvector of h{out). Even if h{out) is hermitian, it has complex 
eigenvalues because we are using a new spectral decomposition, which is only 
possible because we are working in a convenient Gel'fand triplet. This fact 
will be the main tool that we will use below. The eigenvalues and their 
squared will be written as: 



fin = ^^n- -In, In > 



(57) 



m^a^outf' + Zn (58) 

compared with eq. (|34D this equation shows that there are also poles in 
the variable fc^ which are complex numbers. 
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3 Time evolution. 



Coming back to the WKB expansion 
der and the time defined in eq. 
Schroedinger equation: 



if we now consider the next or- 
the function x(a, ip) must satisfy the 



(59) 



where h{ri) is hamiltonian h written as a function of 1]. Even if this hamil- 
tonian is time depending we can consider that for scales a > aout there 
is not particles creation and therefore we have an invariant adiabatic vac- 
uum \0,out > and a definitive pole structure for the S-matrix so for 
a{ri) > ttout expansion ( p^ will always have the same structure. Thus the 
time evolution of x be: 



Xiv) = exp[-i J h{ri)dri]xiO) 



(60) 



From this equation we can obtain some conclusions: 

i.- In particular the time evolution of the right eigenvector \n+ > reads: 



\n{ri)+ >= exp[—i J r2„(r7)(i?7]|ra(0) 



+ > 



(61) 



because even if the pole structure remains fixed the poles move as can be seen 
in the example of ref. eq. (3.3). So, from eq. (|57|) we can see that, if 
there are some 7„ > 0, some of these eigenvectors have a dumped evolution. 
Therefore, these eigenvectors correspond to decaying states. Thus our for- 
malism naturally yields decaying states that vanish towards, the direction of 
time, that we can call the future. 

ii.- Using eq. ( ^2]) we can expand any function \ip >E 0_as: 



(62) 



\ip >= \n+ >< n — \(p > + / >< lu — \(p > duj 

JO 

then its time evolution will be: 

\(p{ri) >= exp[— z / ^ln{ri)dri]\n+ >< n — \ip > + 

n 

roc r 

+ J exp[— 2 J ilk{v)dv]\^+ X UJ — \(p > du 



(63) 
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where all terms in the summatory, such that 7^ 7^ 0, have a decaying evolu- 
tion, while the rest of the terms and the integral have an oscillatory behavior. 
So in the time evolution of (almost) any state we have a decaying term that 
vanishes towards the future. 

iii.- In this way the asymmetry introduced in hypothesis H2 produces an 
effective time asymmetry, because it allows to define a future time direction. 
Moreover, it can be proved that, if in eq. ( |60|) the evolution operator is 
considered as an operator from space 0_ to space 0_, namely if we restrict 



the dynamics to space 0_, eq. (|60|) is only defined for rj > 0. Therefore 
the evolution operator cannot be inverted and so, it is really an irreversible 
operator {see [^, 0], [0). 



iv.- Let us consider the case of mixed states. For a mixed state p G $_ 
we can generalize the spectral decomposition (|62D to obtain: 



n+ >< n + \ + 'S2 / Pnuj\'n+ >< uj + \duj+ 
„ JO 



n,m n 



roc roo roo 

+ / y2 P'-'n\uJ+ >< n + \duj + / / p^^'\uj+ >< UJ + \dujduj' (64) 

Jo n JO Jo 

Repeating the computation of the pure state case, we can compute the time 
evolution of state p{ri). Since either 7^ = or 7„ > there will be oscillating 
terms and dumped ones then we obtain: 

P{v) = P*{v) + J ldv]pi{v) (65) 

where the first term, of the r.h.s. is an oscillatory term and the second a 
decaying term, where we have written a first factor corresponding to the 
slowest dumping factor, namely 7 is the smallest of the non-zero 7„. When 
?7 — 00 we have: 

p{v) P*{v) (66) 

P^It]) is a thermodynamical equilibrium state. In fact, since in its evolu- 
tion there are not dumping factors, it behaves like an ordinary stable quan- 
tum state and its entropy is time-constant, namely the one that corresponds 
to thermodynamical equilibrium (below we will normalize this constant to 
zero). It is logically a non stationary oscillatory equilibrium state, because, 
even if it is in thermic equilibrium, the field cannot go to dynamical equi- 
librium since, in our simple model, there are no interaction terms among 
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the field components. If tliese terms would be present, new dumping factors 
would also be present and the final equilibrium would be a stationary state 

13- 

This is the essence of our formalism. Below we will see the results that 
we can obtain if we follow this road. 



4 Decoherence and correlations. 



In reference using our formalism, it is proved that, if the S-matrix has an 
infinite number of complex poles, we have decoherence and that, in unstable 
states, configuration and momentum are correlated, in such a way that the 
universe ends in a classical phase. In this demonstration hypothesis H3 plays 
an essential role. In paper it was not proved that, in general, the S- 
matrix, relevant for our problem, has an infinite set of complex poles, but 
that set was computed in one example, while other examples were proposed. 

Here we will not repeat the demonstration of ref. |2^, but we will com- 
plete this paper observing that using a potential, with a bounded support, as 
in the present paper, the existence of an infinite set of poles is a consequence 
of ref. [|7| and [ID. 



In fact: A massive scalar field, conformally coupled, in metric (^Tj) satisfies 
Klein-Gordon equation: 

(V^V^ + + = 0, R = 6a-^0 (67) 

This equation leads, by variable separation to: 

th = o f(ri) exr)(±ikx) (68) 

(271)1 a{r])^" ^ ^ ^ 

and f{ri) satisfies a generalized oscillator equation with time-depending fre- 
quency: 

/" + u;^{ri)f = 0, ....LU = [a\ri)m'^ + k^]^ (69) 

On the other hand, in ordinary quantum mechanics, the stationary solution 
u{x) of a massive particle in a potential W{x) satisfies: 

m" (x) + uj'^{x)u{x) = 0, uj{x) = [2m{E - W{x))]^ (70) 
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So both phenomena can be mathematically related according to the analogy: 



a^{ri)m^ ^ ~2mW{x), k'^ ^ 2mE (71) 

More details about this analogy can be obtained from paper [p7[] . 



Now from reference page. 218, we know that the S-matrix of a cut-off 
potential W{x), namely a potential with a bounded support, has an infinite 
number of complex poles. Our potential is ~a^(?7) which is practically a 
constant for a ^ ai, since e = 0, and exactly a constant if we consider the 
real value e = so, subtracting this constant final value, we can say that it 
is a cut-off potential, with a S-matrix endowed with an infinite set of poles. 

Almost all potentials used in the literature of Quantum Field theory in 
Curved Space-time are very well behaved in the infinities and can be 
mimicked by this, bounded support, kind of potentials (this is not the case for 
some QC potentials, that we will discuss in the conclusions). So the existence 
of an infinite set of poles seems quite a general feature of the theory. Thus, 
using the equation of paper ||2^, it can be proved that our formalism leads 
to decoherence, to correlations, and to the outcome of a classical universe, in 
a natural way. 

Finally in particular subsystems of the universe the S-matrix has poles if 
unstable quantum states exist in the subsystem [^. Of course, these poles 
will also appear in any complete S-matrix of the universe. 



5 Entropy 

Let p{7]) be the density matrix, of the universe or one of its subsystems, 
for a physical admissible state (p G $_) and let p*(?7) be the corresponding 
thermodynamical equilibrium matrix. In the universe these matrices are 
related by eq. (|65|) . The hamiltonian of the subsystem is necessarily a term 
of the general hamiltonian of eq. (P5|) and its S-matrix must have poles if 



the subsystem is not trivial (see papers and [0]). So we can repeat all 
what we have said for the universe for the case of the subsystem and we can 
also choose a t-asymmetric regular space state for the subsystem. We must 
take care that the dumping or future direction of the subsystem coincides 
with the dumping or future direction of the universe for consistency. I. e. 
the local and global arrows of time must coincide. Then we will also find eq. 
(pSf) for the subsystem. The only difference would be that, if the subsystem 
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hamiltonian h is no time dependent (with respect e. g. to the proper time 
t) in the integral of eq. (|65D must be subtituted by the usual product ht. 

Then we can define the conditional entropy S = S[p{ri)\p^{ri)], of states 
p{ri) with respect to state p*{ri) both for the universe or the subsystem: 

S = S[piv)\pM] = -tr{pir^)\og[p;\v)piv)]} (72) 

such that S[p^,{r])\p^{ri)] = 0, namely the entropy vanishes at equilibrium. 
We well, use the corresponding classical definition, since really we are 

interested in the classical phase of the universe, but, in order to use only one 
notation, we will use the quantum formulae. 
From eq. (|65D we see that: 

trp{r]) = trp^{ri) = 1 ^ trpi{ri) = (73) 

i.e.: if the states p{ri) and p*{ri) are normalized as it should be, pi(?7) has a 
vanishing trace, thus pi{ri) is not a state but the coefficient of a correction 
of the equilibrium state to obtain state p{ri). The vanishing of the trace of 
pi is directly proved in papers ||3^,[^, [^, using our formalism. Now if we 



expand the logarithm in eq. (0), and use eqs. ( [73D we obtain: 

S = S[piv)\pM] = -exp[-i / ^dr^]tr[p:\v)plir^)] + ... (74) 



where the dots symbolize higher order terms 
This entropy has the property: 

\imS[piv)\pM] = (75) 

namely the entropy evolves towards its null equilibrium value. This is so 
because the prefactor in eq. (|7^ dominates any other time variation, since 
p*(?7) is usually oscillatory (namely it will be oscillatory in the case a, but 
not in the case b, see below) and pi{ri) has oscillatory terms and dumping 
factors that vanish faster than the dominant decaying factor. 

In eq. ( |72D we have two matrices: p{ri) and p*{ri), then we have also two 
possibilities. Either both matrixes have the same kind of evolution or they 
have different ones. The first is the general case, but the second case appears 
when p*(?7) follows a different evolution due to, e.g., an external agency. Let 
us consider the two cases: 
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a.- Both matrices follow the same evolution law. Namely, if we have a 
time variable hamiltonian /i(?7), as in the case of the universe, the evolution 
will be: 

P{v) = exp[-z / l{r])dr]]p{r] ), p,{r]) = exp[-z / l{r])dr]]p^{T] ) (76) 

Jr^' Jrj' 

where l{ri) is the corresponding Liouville operator, i. e., l{ri)p = h{i])p — 
ph{ri). 

In subsystems of the universe where the hamiltonian h is not a proper- 
time variable( e. g. in subsystem which does not expand or contract due to 
an external agency), we would have: 

p{t) = exp[-il{t - t')]p{t') p,{t) = exp[-il{t - t'))]p,{t') (77) 

where / is the Liouville operator corresponding to hamiltonian h and t is 
the proper time. In this case the equilibrium matrix p^, (t) also evolves in the 
same way than matrix p{t). Now, from what we have said in Sec. 3, point iii, 
eq. (|76D is only valid if rj > rj'. Analogously, if the subsystem S-matrix has 
poles and we have also chosen for the subsystem a similar admissible function 
space eq.([77D is only valid if t — t' > ot t > t' . Therefore the last two 
evolutions are irreversible. Now, from ref. (using the classical-quantum 
analogy since we are in the classical phase, the "exp " operator of the last two 
equations will be a Frobenius-Perron operator, and we can use the classical 
definition of conditional entropy) we know that: 

S[pmpM] > S[p{r^')\p,iv% ...S[p{t)M)] > SW)M')] (78) 

respectively. It would be = if the evolution operators would be reversible, 
but they are not (consider also eq.( ^4|)). Then these entropies are really 
monotonically growing. Therefore we have proved the second law of thermo- 
dynamics, for the whole universe or for any non trivial subsystem. So our 
formalism yields this fundamental law naturally (compare with the much 
more complicated coarse graining method of paper ||45[|). 

The demonstration is based in the fact that (in both cases) p is an admis- 
sible state (like the ink drop spreading in the glass of water), so p G $_ and 
7n > 0. If we would have taken p G it would be 7„ < and the entropy 
would decrease showing that: 

i.- In the case of the subsystem. $_|_ is the space of non admissible solution 
(the ink drop contracting spontaneously). In fact, in this case the arrow of 
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time is the one of the universe and not the one of the subsystem, and in the 
subsystem we will see a decaying of the entropy, showing that the state is 
not physically admissible. 

ii.- In the case of the universe. Going from $_ to $+ we have simply 
changed our convention (since all possible arrows of time are embodied in 
the universe), calling the "future" the direction of decreasing entropy. 

b.- Let us now consider the important case of a subsystem of the universe: 
the matter and radiation within an expanding universe ( i. e. we do not take 
into account the entropy of the gravitational field, and we consider that this 
filed as an external agency that expands the space, where the matter and 
the radiation are located). Then the conditional entropy is not necessarily 
monotonically increasing, at least for short times. In fact, we cannot use eq. 
( |78D since p^irf) does not satisfy an equation like (|76|2), because its evolution 
is fixed by the external universe expansion. This is completely logical since 
S[p{'n)\P*iji)\ is just the matter radiation entropy (with no gravitational field 
entropy contribution) in an expanding (or contracting) universe with a equi- 
librium state p*{i]), that varies in an independent way and it is well known 
that matter (let us say a gas) can have decreasing entropy into a variable 
geometry (let us say a box with moving walls). A phenomenological study of 



the problem can be found in ref. In this case what we have called 

up to now S is just the entropy gap AS with respect to a variable maxi- 
mal possible entropy S'max- The actual entropy, which grows monotonically, 
is = Sact = "S'max + AS . But AS has not this property. Furthermore, the 
diminishing of AS for short times is welcome, as we will see in the next two 
sections. 

Finally let us consider the origin of the "miracle" that allows us to define 
a growing entropy with no coarse-graining. The miracle is produced by the 
generalized spectral decompositions (p^ ) or (0) which yield the evolution 
(|63D with only dumping (and oscillatory) factors, which, in turn, are only 
possible if we introduce an asymmetry, like the one of hypothesis H2. Namely 
if we restrict the regular state space in a t-asymmetric way we obtain also 
the asymmetric generalized states of expansions and eqs. (^) and (|5^. Is 
this fact so strange? Experimentally we know that we only have a finite set 
of physical measurements to fix a quantum state, e. g. we know the value 
of the wave function in a finite set of point of the configuration space Q. 
In the limit we will have a set of discrete data that allows us to define by 
interpolation e. g. a polynomial belonging to the set of polynomials V of 
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the introduction. But in order to work we need the whole wave function 
complete superspace and therefore we must complete the space expanded by 
the V. In the ordinary theory we complete the space expanded by the V with 
the topology of the norm of Ti or the nuclear one of S. Up to this point we 
have a time-symmetric theory. If we want to introduce time asymmetry we 
coarse-grain the system, to obtain a new complete space of relevant states. 
In the new theory we complete using Hardy-Schwartz function (i. e. we 
complete with the nuclear topology of So we only use one step, instead 
of two, and obtain the same physical results. But, really both methods have 
essentially the same physical base: we have only a finite (or discrete) amount 
of information and this information must be worked out, somehow, to obtain 
the complete space of regular states of the theory. So, as expected, there is 
no miracle. 



6 The entropy gap. 

In this section we study the universe entropy gap AS" = — 5'max, following 
a qualitative idea of Paul Davies |H6|. We will complete this idea actually 



computing the entropy gap after decoupling time. Therefore we change our 
model, it will be still homogeneous and isotropic, with metric ( ^1]), but 
obviously the particle production is finished, so we will consider that we are 
simply in a flat geometry, matter dominated, universe. 

It is well known that the isotropic and homogeneous expansion of the 



universe is a reversible process with constant entropy In this case the 
matter and the radiation of the universe are in a thermic equilibrium state 
p*(t) at any time t. As the radiation is the only important component, from 
the thermodynamical point of view, we can choose p*(t) as a black-body 
radiation state [|^, i. e. p^{t) will be a diagonal matrix with main diagonal: 

p^iuj) = ZT-''^- (79) 
er — 1 

where T is the temperature, uj the energy, and Z a normalization constant 
O, eqs. (60.4) and (60.10)). The total entropy is: 

1 a 

S = yfxKT^ (80) 
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(||50||, eq. (60.13)) where a is the Stefan-Boltzmann constant and V a com- 
moving volume. 

Let us consider our isotropic and homogeneous model of universe with 
scale a. Any commoving volume evolves as V ~ a^, and, since from the 
conservation of the energy-momentum tensor and radiation state equation, 
we know that T ~ a~^, we can verify that S = const. Thus the irreversible 
nature of the universe evolution is not produced by the universe expansion, 
even if p*(t) has a slow time variation. 

Therefore, after decoupling time, the main process that has an irreversible 
nature is the burning of unstable H in the stars (which produces He and, 
after a chain of nuclear reactions, Fe). This unstable state produces poles 
in the corresponding S-matrix and a nuclear reaction process, with mean 
life-time tiyR = 7^^ . Therefore, using eq. (|65D, and considering that 7 
is constant (under proper time variation), since it corresponds to a local 
process considered in the paragraph "a" of the last section, (or simply on 
phenomenological grounds) we can then say that the state of the universe, 
at time t, is: 

p(t) =p.(t)+pie"^* + 0[e-^*] (81) 

where pi is a certain phenomenological coefficient, which is constant in time 
since all the time variation of nuclear reactions is embodied in the exponential 
law e~^^. Also on phenomenological grounds, we can foresee that pi must 
peak strongly around Ui, the characteristic energy of the nuclear process. All 
these reasonable phenomenological facts can also be theoretically explained 
in different ways, e. g.: Eq. (E^ can be computed with the theory of paper 



36] . In reference P0[ it is explicitly proved that pi peaks strongly at the 



energy Ui. So using eq. ([72|) we can compute the entropy gap: 

AS = -tr[plog(p;V)] (82) 

Using now eq. (0), and considering only times t ^ tNR = 7~^ we can 
expand the logarithm, as in eq. (|7^, to obtain: 

AS ^ -e-Hr (p;V?) (83) 

where we have used eq. (|73D . We now introduce the equilibrium state (^) 
for ci; 3> T . Then: 

A5 ^ -Z-^Th-nr{e^pl) (84) 



25 



where is a diagonal matrix with this function as diagonal. But as pi is 
peaked around uji we arrive to a final formula for the entropy gap: 



AS ^ -CT^e-^'e^ (85) 

where C is a positive constant. 

Let us now compute the time evolution of the entropy gap. We have 

2 

computed AS* for times larger than decoupling time and therefore, as a ~ 
and T ~ a~^, we have: 



T = To (^^ j (86) 
where to is the age of the universe and To the present temperature. Then: 

A5^-Cie-^*rtt(^)' (87) 

where Ci is a positive constant. Drawing the corresponding curve it can 
be seen that AS* has a maximum at t = ten ^ minimum at t = tcr2- Let 
us compute these critical times. The time derivative of the entropy reads: 



A5? 



-7 - 2t" 



2 uji ft, 



3 toTo V t 



'^0 



AS (88) 



This equation shows two antagonic effects. The universe expansion effect is 
embodied in the second and third terms in the square brackets, being an 
external agency to the matter-radiation system such that, if we neglect the 
second term, it tries to increase the entropy gap and, therefore, to take the 
system away from equilibrium (as we will see the second term is practically 
negligible). On the other hand, the nuclear reactions embodied in the 7- 
term, try to convey the matter-radiation system towards equilibrium. These 
effects become equal at the critical times tcr such that: 



to 2uJi f t 



7to + 2^ = 5^(^)" (89) 

For almost any reasonable numerical values this equation has two positive 
roots: ten <^ to ^cr2- Precisely: 
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i.- For the first root we can neglect the 7)f:o-term and we obtain: 



(this quantity, with minus sign, gives the third unphysical root). 

ii.- For the second root we can neglect the 2(to/^cr)~term, and we find: 

Let us make now some numerical estimates. We must choose numerical 
values for four parameters: Ui = Tj^r, t^R = 7~^, ^o; and Tq. 



Tj^R and Inr can be chosen between the following values ||51[] : 

Tnj^ = 10^.to..lOl°ir (92) 

t^R = 10^.. to.. 10^. years 
while for to and Tq we can take: 

to = 1.5 X lO^^.years (93) 

To = 3^K 

In order to obtain a reasonable result we choose the lower bounds for Tjvr 
and tjvij and for t^n we obtain : 

ten ~ 1-5 X lO^.years (94) 

So ten is smaller than the decoupling time and it should not be considered 
since the physical processes before this time are different than those we have 
used in our model. Also, we must only consider times t > tNR = 7~^, in 
order to use eq. (^31). 



For tcr2 we obtain: 



tcr, ^ 10% (95) 
From eqs. (|9^ and (0) we can see that really tcri ^ to tcr2- Thus: 



-From tNR to tcr2 the expansion of the universe produces a decreasing of 
entropy gap, according to Paul Davies prediction [0. Also, it probably pro- 
duces a growing order, and therefore the creation of structures like clusters, 
galaxies and stars 



27 



-After tcr2 we have a growing of entropy, a decreasing order and a spread- 
ing of the structures: stars energy is spread in the universe, which ends in a 



thermic equihbrium [Q. In fact, when t — oo the entropy gap vanishes (se 
eq. (0)) and the universe reaches a thermic equihbrium final state. 

tcr2 ^ lO'^to is the frontier between the two periods. Is the order of 
magnitude of tcr2 ^ realistic one? In fact it is, since lO^to ~ 1.5 x lO^'^years 



after the big-bang all the stars will exhaust their fuel ||5^ , so the border 
between the two periods should have this order of magnitude. Further, it 
should also be smaller than this number. This is precisely the result of our 
calculations (see also letter fETH). 



So we are at the edge of a correct physical prediction, even if our model 
is extremely naive and simplified: an homogeneous universe and besides we 
have neglected the higher order terms in eq. ( pT]) which perhaps may be 
important for finite times. Besides in the real universe nuclear reactions take 
place within the stars, that can only be properly considered in a inhomo- 
geneous geometry. Nevertheless, this rough numerical estimate shows the 
consistency of all the theory. Furthermore the decreasing of the entropy gap, 
in the period tNR < t < tcr2, will be crucial in the next section. 



7 The branch system. 

The set of irreversible processes within the universe, each one beginning in 



an unstable non-equilibrium, state can be considered a branch system |46 



2^, . Namely, every one of these processes begin in a non-equilibrium state, 
such that, this state was produced by a previous process of the set. E. 
g.: Gibbs ink drop (initial unstable state) spreading in a glass of water ( 
irreversible process) is only probable (since the probability to create an ink 
drop by fluctuations is extremely small) if there was first an ink factory, which 
extracted the necessary energy from an oven, where coal (initial unstable 
state) was burnt (branched irreversible process); in turn coal was created 
with energy coming from the sun, where H (initial unstable state) is burnt 
into He (branched irreversible process); finally H was created using energy 
obtained from the unstable initial state of the universe (the absolute initial 
state of the branch system). Therefore, using this hierarchical chain, all the 
irreversible processes are related to the cosmological initial condition, the 
only one that must be explained. Let us observe that: 
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i. - The branch system defines its own arrow of time the branch arrow 
of time (BAT), as the direction that goes from the unstable initial state of 
every member of the system towards equilibrium. Probably the BAT is the 
most useful of all arrows of time, since it is present in any irreversible local 
process. 

ii. - Once we have the branch system the irreversible evolution of each 
system is easy to explaining, since once we have understood the origin of the 
initial unstable state of each irreversible process within the universe (even if 
we have not yet discussed the origin of the initial state of the whole universe) 
it is not difficult to obtain Lyapunov variables (or irreversible evolution equa- 
tions), if we consider, e. g. that the subsystems where these processes take 
place are not isolated. If it is so, forces of stochastic nature penetrate from 
the exterior of each subsystem and, it is well known, that if we add stochas- 
tic terms to time-symmetric evolution equation, we obtain time-asymmetric 
ones, yielding Lyapunov variables e. g. a non- decreasing entropy We 
can also consider that each subsystem has an enormous amount of infor- 
mation and we are able to measure, compute, and control a part of this 
information, that we will call relevant. If we neglect the rest of the informa- 
tion, the irrelevant one, we can obtain also irreversible evolution equations 



and Lyapunov variables [@,|]55[. These two procedures can be considered 
within the coarse-graining usual formalism. 

iii.-But of course, if we follow the ideas of this paper, we will use more 
refined mathematical tools, and in each of the subsystem, introduce a model 
similar to the one we have used for the whole universe, as we have done in 
section 5, introducing the hypothesis H2 in each subsystem. (It has already 
been done in papers 0], [^, and the same results are obtained, i. 

e. irreversible evolution equation, Lyapunov variables, etc.) Then we see 
that entropy grows in each subsystem provided the state of the subsystem 
would be chosen among the physically admissible states of space Then 
each subsystem of the branch system begins in an unstable, low entropy, state 
and evolves towards thermal equilibrium. The physical non-admissible states 
(those of space $+) correspond to theoretical evolutions that would only exist 
before the instant of creation of the subsystem (the instant when we put the 
ink in the glass of water ) evolving with decreasing entropy, towards that 
instant (namely the ink drop contracting spontaneously) . These evolutions 
simply do not exist in nature because, before the instant of its creation, 
the subsystem really does not exist as such. Before that instant a different 
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subsystem exits with different evolution laws (the ink factory that creates the 
ink drop). Therefore all the scenario turns out to be realistic and satisfactory. 

iv. - So only one problem is left: Why did the universe begin in an unstable 
low-entropy state? Let us first observe that really we are referring not to 
the "whole" universe (with its gravitational field) but only to the matter- 
radiation subsystem of the universe. 

In the no-time version of the introduction we have postulated Hi and 
H2. Using these hypotheses we have reconstructed time and demonstrated, 
in the sections above, that the universe expansion creates, in its matter- 
radiation subsystem, an entropy gap AS" that takes it out of equilibrium, not 
only at t = 0, but in a long period of its history, since the actual entropy is 
Sact < 5'max- Wc havc also demonstrated that the matter-radiation subsystem 
of the universe evolves to a final state of thermic equilibrium since AS 0, 
when t — s> 00 (cf. eq. ( |75p). So the answer to the only problem left is 
essentially hypothesis H2. Certainly, someone will think that we have solved 
a problem by postulating an axiom, and this is not a very exiting result. But 
if the axiom yields the solutions of many problems, and this is the case of 
H2, the axiom must be welcome. After all, this is the role of axioms. 

v. - Finally we can ask ourselves if, in the perspective of the branch system 
idea, H2 is a natural hypothesis. H2 says that p G $- G 5 ® 5. So first it 
is postulated ,that p G iS ® 5, and therefore p is a smooth function, with 
infinite derivatives, and well behaved in the infinities of the configuration 
space. This part of H2 seems quite natural. Certainly much more natural 
than the two other alternative possibilities: 

a. - p G £, in which case p can be, e. g., a square integrable function of 
7i ® ?i, where in a set of points, the function can take non continuous an 
arbitrary values. What is the physical meaning of this discontinuity? 

b. - p E <S) , namely a distribution, e. g. a delta function, certainly 
a quite unnatural state. 

So the first part of H2 is natural. The second part is to ask why p would 
be endowed with a natural asymmetry, the one of Is it too much to 
ask? Let us study this question according to the branch system idea and our 
formalism. There will be no branch system only if the universe (and now 
we are referring to the whole universe with the gravitational field included) 
would begin in an equilibrium state p^,, since in this case it will always remain 
in equilibrium. Now, from eqs. (|6^ , (|6^), and (|66|) we see that, in this case 
p* G so p* would be a distribution [Q, something like a delta function. 
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and we have just considered this choice as unnatural. On the contrary p G $_ 
is a much more natural state. Any state of $_ will produce a branch system, 
since any state of $_ yields eq. So we can, at least, conclude that 

H2 is the requirement that the state of the universe id a natural and an 
asymmetric one. H2 is also intimately related with the branch system idea 
and in consequence it is also related with the fact that really our universe is 
a branch system. H2 is just the transcription of these physical facts. 

8 Coordination of the arrows of time. 

In this section we will only consider the coordination of the arrows of time 
related with our model, namely: 

• The branch arrow of time (BAT), the arrow that goes from the unsta- 
ble initial state of every process of the universe branch system to its 
equilibrium final state. As we have seen in the last section this arrow 
is a direct consequence of the asymmetry introduced by H2. 

• The thermodynamic arrow of time (TAT), that points to the direction 
of the growing of the universe entropy S. 

• The cosmological arrow of time (CAT), that points to the direction of 
the growing of the universe scale a. 

Of course, all these arrows are related with time and therefore they must 
only be considered in the classical (or semiclassical) period where the time 
Tj, given by eq. (p9D, is well defined. In the timeless quantum period we only 
have the asymmetry defined by II2. 

Then: eq. (|7Sp, for non-expanding or contacting subsystems within the 
universe, which is a consequence of H2, shows that BAT=TAT, and that t, 
or more generally t], grows in the same direction than S. 

The relation between BAT(=TAT) and CAT is given by the ± sign in eqs. 
(|30|) or (0). Then these two arrows of time, a priori, are not coordinated in 
our model. But in the classical period we have just one classical universe and 
therefore the ± sign and the constant C are fixed, so in the classical period 
we have just one sign: either + or - (see App. 12). Therefore, once the sign 
is fixed, a clear relation appears between BAT and CAT: 



31 



• If the model is an expanding one (and we choose the sign +) we will 
have BAT=CAT, at least in the final evolution where eq. (0) is valid 
(if we make the unusual choice of the - sign we are just changing the 
conventional direction of time 77, with no physical consequences). 

• If the model is an expanding-contracting one (and we choose the sing 
+) we will have BAT=CAT in the expanding period and BATt^CAT 
in the contracting period. BAT=CAT is, in fact, the definition of the 
expanding period and BAT 7^ CAT is the definition of the contracting 
one. But BAT=TAT does not change when we go from the expanding 
to the contracting period or vice versa, since the choice of $_ (or $+) 
is made once and for all. 

So the study of the correlation of the arrows of time is completed, and 
almost trivial, because we have H2 that defines BAT. (See also [p7|] ) 



9 Other results. 



The main results related with quantum cosmology are stated in the above 
sections. But we must comment that using the present formalism all the 
relevant results of irreversible statistical mechanics can also be obtained, e. 
g. all the results of the book [^], as it is proved in ref. ||3^, because the main 
n projector of the quoted book can be defined using Gel'fand triplets. Also, 
in some simple cases, we can go from the quantum models to the classical ones 
where we find the same philosophy, in classical cases. Chaotic models 
like Baker's transformation and Renyi's maps, are also treated with the same 
method, with good results |60|. Other interesting results are contained in 

and 



papers |T^, [0,^9[, , and [^. So what we have explained is just the 
quantum cosmological chapter of a general method to deal with irreversible 
processes. 



10 Conclusions. 

Let us summarize our main conclusions. 

i.- All our scheme is based in the existence of a physically admissible state 
superspace $_ and of a physically forbidden state superspace Thus, the 
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time inversion that goes from to $+ is also forbidden. Namely, no Maxwell 
demon can change the direction of all the velocities of the universe. This is, 
of course a practically impossible task. Is it also theoretically impossible? 
In fact it is, even if Maxwell demon would change all these directions (while 
we are sleeping) we will not notice the change (when we wake up), since all 
arrows of time would be changed and we would not have any extra arrow to 
verify the change. Thus this global-demon task is conceptually impossible. 

ii.- What we have presented is not a mathematical theorem, but a model 
that can be generalized in many ways. These more general models will have 
a similar behavior, than the present one, if two essential features are present: 
the existence, at the quantum gravity level, of a geometry of maximal relative 
probability which allows us to construct "out" states for the fields, and a S- 
matrix with infinite complex poles. The first requirement seems natural for 
any realistic model of universe. On the other hand we have restricted the class 
of possible potentials in order to be able to prove that the corresponding S- 
matrix has infinite complex poles. But several QC potentials do not belong 
to this class, because they have a bad behavior at infinity. Nevertheless, 
usually, they also have an infinite set of poles, as it can be proved case by 
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So the two basic features seem usual enough to consider that our 
model is a good sample of the general behavior of the universe. Then we can 
say that: 

iii. - If we introduce an adequate regular state space (or an adequate topol- 
ogy) it seems that all the known results of statistical irreversible physics can 
be reobtained. It must be emphasized that we are not adding a new object to 
the theory, since a regular state space (or the corresponding topology) must 
be defined anyhow. We are just choosing the most convenient one. Let us 
repeat the general relativity example: The space-time has a metric, we can 
choose a flat space-time metric or a curved one. In the second case we ex- 
plain gravity in the best possible way. We add nothing, we just choose the 
best mathematical structure. The same thing happens in the present case. 
If we choose the usual regular state space S we are forced to make a coarse 
graining (and there is nothing experimentally wrong with coarse graining, 
as there in nothing experimentally wrong with post-newtonian theories, but 
both are "non economical" formalisms). If we choose the new regular state 
space 0_ we make two steps in one, so we have a conceptual advantage. 

iv. - Precisely, because the new formalism is conceptually clearer we can 
see that time-asymmetry is just a kind of spontaneous symmetry breaking. 
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v.- Most probably the old and new formalism will always yield the same 
physical results, because they are both based in the same physical base: the 
limited amount of information must be somehow worked out to obtain a com- 
plete theory. Therefore, most likely, they are as experimentally equivalent 
as general relativity and post-newtonian gravity with an infinite number of 
terms. 

So, even if we have not found any new or spectacular result, we think 
that the introduced formalism presents a quite coherent picture of the real 
time-asymmetric universe and shows us how time-asymmetry forces us to 
choose a Gel'fand triplet as the mathematical structure of the theory. 

11 Appendix. Time-asymmetry of the regu- 
lar states superspace. 

It is a common lore among physicists that, since the fundamental evolution 
equations are time-symmetric and the universe is time asymmetric, therefore 
the state of the universe must be the responsible for time asymmetry. Namely 
time-asymmetry is created as in the case of spontaneous symmetry breaking. 
The problem is to mathematically precise this idea. That is what we have 
done in this paper, by showing that the asymmetry can be hidden in the 
choice of an asymmetric regular state superspace. In fact, beginning with 
the one symmetric regular state space S of the primitive theory, we have 
introduced the asymmetry related with the lower and upper halfplanes, of 
the complex plane, and we have defined asymmetric regular state spaces 
</)„ and </)+, and we have chosen one of them, say (p-, creating a symmetry 
breaking. 

But, since this is perhaps the most important point, we will try to explain 
this fact in the most didactical way. First let us consider that a > 0, and 
let us draw the potential V{r]) = V^[a(?7)] in the right side of the horizontal 
7] axis. Let us take C such that rj = would correspond to a = in eq. 
(|3TD, to make the origin of a to coincide with the one of t]. As, in principle, 
7] can be positive or negative we will complete the previous figure with a 
symmetrical curve in the left side of the vertical axis. Really as we cannot go 
from negative t] to positive ones we can add an infinite potential barrier at 
1] = 0. Now we have a more or less typical potential figure for a spontaneous 
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symmetry breaking. In our case the minima are unimportant (even if, most 
likely, the equilibrium states will correspond to the minima in almost every 
model) but since we have an infinite potential barrier at = it is clear that 
the universe can either be in the right or the left side of the vertical axis. 
Choose either one or the other possibility. Then, you have chosen the + or 
the - sign in eqs. (|30D or (^Tj) (since we have chosen a = when 77 = and 
it is always a > 0). Let us consider, for a moment the expanding case: if 
you choose + then BAT=CAT and you have chosen the convention of all the 
formulae of this paper, namely you have chosen a H2 with a $_ space as in 
this paper. If you choose the - sign, you just change the convention used in 
this paper, instead of rj you will have —r] everywhere and in H2 you will have 
$+. Of course, physically nothing is changed since you have just changed a 
convention. In the expanding-contracting case we can repeat the reasoning 
for the expanding phase only and everything turns out to be the same. 

So we are clearly in a case of a symmetry breaking, produced by the 
state space and not by the equations, where we have conventionally chosen 
between two conventional different possibilities, either the right side of the 
figure or the left side, namely between $_ and $+. After this choice the 
time-asymmetry appears but, of course, the choice itself is irrelevant, as in 
every spontaneous symmetry breaking. So really there are two ± choices in 
the theory: the one that corresponds to the choice between $_ and and 
the one in eqs. (pOD and (|3TD. We will explain more about the second choice 
in the next appendix. 

12 Appendix. The asymmetry of Q space. 

Why do we have + and not - in eqs. (0) or ( |3TD ? Perhaps it is interesting to 
relate the present formalisms with the final speculations of the second paper 

If $_ is chosen (namely the first ± choice is made)and C > : Are 
the + and -, of eq. (^Tj) conventionally identical? Is the choice between 
them just conventional or is there something substantially different between 
them? Of course the difference between + and - is just impossible to find in 
the quantum gravity timeless period. But the arena of time asymmetry is the 
classical period where the ± sign and the constant C are fixed (and as in this 
case we have taken C > 0, as we have said). Can we find a difference between 
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+ and - in this period? Our strategy is to maintain the symmetry the WDW 
equation and create the time-asymmetry using the states. In doing so we 
can use all the features of the states and the configuration space Q. And it 
turns out that Q is not so symmetric since it has boundaries that gives Q 
some peculiar structure. Precisely a > in our very simple model. So once 
the choice of the space of state is made, using all the well known philosophy 
of symmetry breaking, as explained in section Appendix 11, + and - become 
obviously different in expanding models since: 

• With the + we would have equilibrium towards a — »• oo, which is 
possible. 

• With the - we would have equilibrium towards a ^ — oo, which is 
impossible since a > 0. 

So we must choose the +. If we would have begun choosing $4. (and 
C < 0) we will be forced to choose the - sign. So the choice of the regular 
state space implies the choice of a sign. In expanding-contacting models 
things are more difficult since we do not have either a— >-ooora— >— cx), so 
theoretically we never reach the equilibrium state. In these models we must 
change the notion of equilibrium, and consider that this state is obtained in 
long but finite times. Then, we must consider only the expanding phase and 
repeat the reasoning of expanding models. 

Therefore, as it was already anticipated in ref. the boundary structure 
or Q creates a difference, and allows us to choose the right sign. 
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